Dynamic self-assembly and patterns in electrostatically driven granular media 
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We show that granular media consisting of metallic microparticles immersed in a poorly conduct- 
ing liquid in strong DC electric field self-assemble a rich variety of novel phases. These phases include 
static precipitate: honeycombs and Wigner crystals; and novel dynamic condensate: toroidal vor- 
tices and pulsating rings. The observed structures are explained by the interplay between charged 
granular gas and electrohydrodynamic convective flows in the liquid. 
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Understanding the unifying principles of self-assembly 
of complex systems such as macro- molecules 0, diblock 
copolymers [^), micro-magnetic systems ensembles 
of charged particles [[| is the key to future advances in 
nanoscience. Large ensembles of small particles display 
fascinating collective behavior when they acquire an elec- 
tric charge and respond to competing long-range electro- 
magnetic and short-range contact forces. Many indus- 
trial technologies face the challenge of assembling and 
separating such single- or multi-component micro and 
nano- size ensembles. The dynamics of conducting mi- 
croparticles in electric field in the air was studied in j|, |) . 
Phase transitions and clustering instability of the electro- 
statically driven granular gas were found. The studies of 
self- assembly of colloidal particles in aqueous solutions 
revealed the importance of self-induced electrohydrody- 
namic (EHD) convective flows on the formation of vari- 
ous precipitate states 

B % E !> 0- O rdered clusters of 
particles vibrated in liquid were studied in Ref. Q . 

In this Letter we report new dynamic phenomena oc- 
curring in granular gas in poorly conducting liquid sub- 
ject to strong electric field (up to 20 kV/cm). We show 
that metallic particles (120 /im diameter Bronze spheres) 
immersed in a toluene-ethanol mixture in DC electric 
field self- assemble into a rich variety of novel phases. 
These phases include static precipitates: honeycombs 
and Wigner crystals; and novel dynamic condensates: 
toroidal vortices and pulsating rings (Figs, [j] and ||). 
The observed phenomena are attributed to interaction 
between particles and EHD flows produced by the action 
of the electric field on ionic charges in the bulk of liq- 
uid. This provides a new mechanism for self-assembly of 
microparticles in non-aqueous solutions. 

To form the electro-cell, granular media consisting of 
3 g (about 0.5 x 10 6 ) mono-dispersed 120 /im bronze 
spheres was placed into a 1.5 mm gap between two hor- 
izontal 12.5x12.5 cm glass plates covered by a transpar- 
ent conducting layer of indium tin-dioxide (the particles 
constitute less than a monolayer coverage on the bot- 
tom plate). An electric field perpendicular to the plates 
was created by a DC high voltage source (0-3 kV) con- 
nected to the inner surface of each plate [|| || . The liquid 




FIG. 1: Snapshots of a) static clusters; b) honeycomb precip- 
itate c) "down" toroidal vortices; d) pulsating rings. 




FIG. 2: Close-up view of (a 
(b) honeycomb precipitate 



Wigner crystal precipitate and 



was introduced into the cell through two teflon micro- 
capillaries connected to opposing side walls of the cell. 
Real time images were acquired using a high speed, up 
to 1000 frames per second, digital camera suspended over 
the transparent glass plates of the cell. Experiments were 
repeated with smaller particles (40 fim Cu spheres) and 
qualitatively similar results were obtained. 
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The principle of the electro-cell is as follows: a parti- 
cle acquires an electric charge when it is in contact with 
the bottom conducting plate. It then experiences a force 
from the electric field between the plates. If the upward 
force induced by the electric field exceeds gravity, the par- 
ticle travels to the upper plate, reverses charge upon con- 
tact, and is repelled down to the bottom plate. This pro- 
cess repeats in a cyclical fashion. In an air-filled cell, the 
particle remains immobile at the bottom plate if the elec- 
tric field E is smaller than a first critical field E\. If the 
field is larger than a second critical field value, E2 > E\, 
the system of particles transforms into a gas-like phase. 
When the field is decreased below E 2 (E 1 < E < E 2 ) lo- 
calized clusters of immobile particles spontaneously nu- 
cleate to form a static precipitate on the bottom plate. 
The clusters exhibit the Ostwald-type ripening jl2[ . The 
phase diagram for the air-filled cell is symmetric with 
respect to the reversal of the electric field direction. 

In this work, we discovered that the situation is re- 
markably different when the cell is filled with a non- 
polar low-viscosity liquid. For the liquid, we typically 
used HPLC grade toluene with various concentrations of 
ethanol (98% w/w anhydrous ethanol + 2% w/w water). 
The conductivity of the liquid changes from 5 x 10 -11 
Ohm -1 m _1 (pure toluene) to 5 x 10 -9 Ohm -1 m _1 
(9% ethanol solute). The phase diagram delineating the 
different particle behavior as a function of the applied 
voltage and concentration of ethanol is shown in Fig. ||. 
For relatively low concentrations of ethanol solute (be- 
low 3%), the qualitative behavior of the liquid- filled cell 
is not very different from that of the air-filled cell: clus- 
tering of immobile particles and coarsening are observed 
between two critical field values E\^ with the clusters be- 
ing qualitatively similar to that of the air cell. A typical 
low-concentration pattern is shown in Fig. |l|a. 

When the ethanol concentration is increased, the be- 
havior of the liquid-filled cell is remarkably different from 
that of the air-filled cell. The phase boundaries are no 
longer symmetric with respect to reversal of the applied 
voltage. Initially, an increase in ethanol concentration 
leads to asymmetric behavior with respect to the direc- 
tion of the electric field. Critical field values, £1,2, are 
larger when the electric field is directed downward (" +" 
on upper plate) and smaller when the field is directed 
upward ("-" on upper plate). This difference grows with 
increasing ethanol concentration. 

The observed asymmetry of the critical fields is appar- 
ently due to an excess negative charge in the bulk of the 
liquid. This charge increases the effective electric field 
acting on the positively-charged particles sitting on the 
bottom plate when a positive voltage is applied to that 
plate and correspondingly decreases the effective electric 
field in the case of an applied negative voltage. From the 
shift in the critical field values one estimates the excess 
negative charge concentration in the bulk. Assuming uni- 
form charge distribution in the bulk one derives the esti- 
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FIG. 3: Phase diagram, U is applied voltage (positive val- 
ues correspond to plus on top plate), C is concentration of 
ethanol. Domain sc denotes static clusters, w-Wigner crystal, 
hc-honeycombs, utv and dtv up/down toroidal vortices. 



mate for the excess charge concentration n ~ esq AU / ecP, 
where the dielectric permeability of the toluene/ethanol 
mixture e w 3, d is the gap between the plates, and 
e is ionic charge. For the solute concentration C=3 % 
and the observed critical voltage shift AU = 150 V, one 
finds n~ — n + ~ 10 16 m -3 . The carrier concentration can 
be independently estimated from transport measurement 
(total current through the cell is / ~ 1 /aA at U = 300 
V, its surface is S = 1.5 x 10 -2 m 2 ). Neglecting concen- 
tration gradients, the migration current can evaluated as 
e 2 DS(n~ + n+)U/dk B T, where D « 10~ 9 m 2 /s is the 
representative value for ion diffusion, and T is temper- 
ature. That estimate also yields n~ + n + ~ 10 16 m -3 . 
Both estimates indicate that the liquid contains almost 
exclusively ions of the same (negative) sign (albeit in a 
very small concentration). These ions (possibly OH~) 
are distributed throughout the cell by the current, while 
positive ions (i^ + ) are absorbed by the electrodes |l3| ]. 

The situation changes dramatically for higher ethanol 
concentrations: increasing the applied voltage leads to 
the formation of two novel immobile phases: honeycomb 
and Wigner crystal precipitates. The Wigner crystal pre- 
cipitate (Fig. ga) is induced for upward direction of the 
applied electric field, whereas a honeycomb precipitate 
is induced upon electric field reversal (Figs. [j]b, gb). A 
return to the original applied voltage transforms the hon- 
eycomb structure back into the Wigner crystal phase and 
vice versa [Q. The representative time scale for these 
transitions is 0.5 sec. For high enough voltage these pat- 
terns coexist with granular gas. However, due to the 
bombardment by the gas- phase particles, the honeycomb 
structure becomes more diffuse. The typical size of the 
elementary honeycomb cell is determined by the gap be- 
tween the plates. Increasing the gap from 1.5 mm to 3 
mm results in doubling of the honeycomb cell size. 

A further increase of ethanol concentration leads to the 
appearance of a novel dynamic phase - condensate (Fig. 
|l|c,d). The dynamic condensate is very different from the 
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FIG. 4: Upper panel: two horizontally rotating down t- 
vortices. Arrows indicate direction of horizontal rotation of 
binary t-vortices. Small arrows show particle trajectories at 
top plate for single vortex. Inset: vertical cross-section of 
vortex flow for down t-vortex. Lower panel: time evolution 
of one rotating asymmetric down t-vortex. 



static precipitate observed for low ethanol concentrations 
and voltages: almost all particles in the condensate are 
engaged in a circular vortex motion in the vertical plane, 
resembling that of particles in Rayleigh-Benard convec- 
tion, see Fig. [| jl(J. The condensate co-exists with the 
dilute granular gas. The direction of rotation is deter- 
mined by the polarity of the applied voltage: particles 
stream towards the center of the condensate near the 
top plate for upward field direction and vice versa. The 
direction of rotation is visible using particle-image ve- 
locimetry. In other words, for upward field direction the 
particles stream down at the center of the condensate 
domain forming downward toroidal vortices (Fig. |l|c). 
Conversely, the particles form upward toroidal vortices 
for downward field direction. 

The evolution of the condensates depends on the elec- 
tric field direction. For downward field, the "up" t- 
vortices become unstable due to the spontaneous for- 
mation of voids in the condensate (Fig. [l|d). These 
voids exhibit complex intermittent dynamics. Medium- 
size clusters typically form only a single void and fasci- 
nating pulsating rings appear (see Fig. ||). In contrast, 
for upward field, large "down" t-vortices remain stable. 
In some cases two or more "down" t-vortices merge into 
one, forming an asymmetric condensate which performs 
composite rotation in the horizontal plane, see Fig. |[ 

The honeycombs, Wigner crystal, and t-vortices are 
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FIG. 5: Time evolution of a pulsating ring. 



likely caused by the same mechanism. It is known 
that electric field generates local EHD convective flow in 
charged liquid in the vicinity of microparticles @, ||, ||. 
We observed that for the downward field direction the 
negatively charged liquid flows upwards around a station- 
ary lying particle, creating an axially symmetric toroidal 
flow field. The upward flow induces a sweeping horizon- 
tal flow on the bottom plate toward the center of the 
precipitate. Such EHD-induced flow can assemble parti- 
cles into a honeycomb structure. The proportionality of 
the honeycombs scale to the gap between plates supports 
the EHD convective mechanism for honeycomb formation 
(compared to Rayleigh-Benard convection). In the case 
of upward field direction, each particle induces a down- 
ward flow of the fluid and, consequently, a sweeping hori- 
zontal flow away from the center of a precipitate, leading 
to mutual repulsion and the formation of a Wigner crys- 
tal. For elevated concentrations of ethanol and higher 
voltages the EHD flow is able to drag bouncing particles 
from the gas phase and assemble them into t-vortex con- 
densates. This mechanism is consistent with the obser- 
vation that particles stream outwards from the center of 
a "down" t-vortex on the bottom plate (as in the Wigner 
crystal) and inward towards the center, of an "up" t- 
vortex (as in honeycombs). Spontaneously- formed voids 
in large, "up" t-vortices are inflated by EHD flow (like 
voids in a honeycomb structure), leading to persistent 
modification of the shape of the pulsating rings. 

The toroidal vortices (t-vortices) attract each other 
and coalesce, resulting in large-scale dynamic conden- 
sates. The representative results illustrating hydrody- 
namic attraction between t-vortices are shown in Fig. ||. 
The separation between two t-vortices, i?, for small dis- 
tances with respect to the cell gap h agrees with the fol- 
lowing power law R ~ (to — t) 1 / 3 , where to = const. This 
behavior yields the following expression for the attraction 
velocity V = dR/dt = —n/R 2 , n = const which strongly 
indicates the hydrodynamic character of the vortex in- 
teraction. Indeed, in our quasi-two-dimensional geome- 
try, the asymptotic hydrodynamic velocity v for inviscous 
flow far away from the vortex behaves as v ~ 1/R. Due to 
incompressibility J vu(z)dz = 0, where z is the vertical 
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FIG. 6: Distance R between two vortices vs time t. Symbols 
show experimental data, solid lines are fits according to Eq. 
(0), dashed lines are R ~ (t -t) 1/3 , for U = -800 V, C = 8% 
w/w (dtv) and U = 1900 V, C = 8% w/w (utv). 



not observed. The main difference here is the presence of 
thick spatially-charged layer in the dielectric liquid. The 
fixed sign excess charge in the bulk results in a different 
kind of interaction between particles and the EHD flows. 

In conclusion, we report the discovery of a series of 
remarkably rich novel vortex phases and the dynamic 
assembly of electric field-driven conductive particles in 
poorly conductive fluids. Our results exemplify the im- 
portance of interactions between the particles and un- 
compensated ionic charges in the bulk of the liquid. Ex- 
tension of our experiments toward self-assembly of submi- 
cron particles may have intriguing practical applications, 
such as lithographic patterning and assembly ||, [l7| , 
etc. We are grateful to A. A. Abrikosov and V.M. Vinokur 
for fruitful discussions. This research was supported by 
the US DOE, grant W-31-109-ENG-38. 



coordinate and v\\ is horizontal hydrodynamic velocity. 
Therefore, in the lowest order, the interaction between 
t-vortices is determined by the square of the hydrody- 
namic velocity vy, i.e. V ~ v 2 ~ — 1/R 2 . In the presence 
of other t-vortices and due to influence of the boundaries 
the interaction law is modified as follows 



dR 
~dt 



R 



K 
R 2 



0(V 2 ) 



(1) 



where Vo = * s mean hydrodynamic velocity cre- 

ated by other vortices and ( = const. Assuming Vo « 
const (which is true when all other vortices are far away 
from the two vortices of the interest) one obtains excel- 
lent fitting for the experimental curve in Fig. ^. 

An order-of-magnitude estimate of EHD velocity can 
be obtained from the balance of viscous (nV 2 v) and elec- 
tric (enE) forces acting on the charged liquid. Here 
n = 5.6 x 10 -4 Pas is the viscosity of toluene and n 
is the ion concentration. Since the velocity varies on a 
characteristic scale of the order of the particle radius a, 
one can replace V 2 with a -2 , and thus v « enEa 2 /n. 
Taking n ~ 10 16 m -3 gives an EHD velocity of ~ 5 
mm/sec at an applied voltage U = 500 V. This estimate 
is in qualitative agreement with the measured times of 
pattern formation and vortex attraction. 

We also carried out experiments with solutions of trib- 
utylamine tetraphenilborate (TBATPB) in toluene. This 
system shows the same phases but with opposite polar- 
ity of the applied voltage. The critical lines in the phase 
diagram of the TBATPB/toluene mixture have slopes 
opposite to that of the ethanol/toluene based phase dia- 
gram (see Fig. |3|). This observation depicting reversal of 
behavior under opposite polarity indicates the presence 
of positive ions in the bulk of the solution. 

There is a significant body of research on the behavior 
of dielectric colloidal particles in aqueous solutions [see 
e.g. Refs. [|| |?], ||, However, dynamic patterns were 
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